Memory erasure in small systems 
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We consider an overdamped nanoparticle in a driven double-well potential as a generic model 
of an erasable one-bit memory. We study in detail the statistics of the heat dissipated during an 
erasure process and show that full erasure may be achieved by dissipating less heat than the Landauer 
bound. We quantify the occurrence of such events and propose a single-particle experiment to verify 
our predictions. Our results show that Landauer's principle has to be generalized at the nanoscale 
to accommodate heat fluctuations. 
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Maxwell's demon is a salient figure of thermodynamics 
Introduced in 1867 to illustrate the statistical nature 
of the second law, the demon is an intelligent microscopic 
being able to monitor individual molecules contained in 
two neighboring chambers initially at the same tempera- 
ture. By opening and closing a small hole in the separat- 
ing wall, the demon collects the faster ('hot') molecules 
in one of the chambers and the slower ('cold') ones in 
the other, thus creating a temperature difference. The 
demon is therefore able to decrease the entropy of the 
system without performing any work, in apparent viola- 
tion of the second law of thermodynamics. The paradox 
was eventually resolved by Bennett who noted that 
during a full thermodynamic cycle, the memory of the 
demon, used to record the coordinates of each molecule, 
has to be reset to its initial state. According to Lan- 
dauer's principle memory erasure necessarily requires 
dissipation of entropy: the cost of erasing one bit of in- 
formation is at least l^Srcmdauer = A: In 2, where k is the 
Boltzmann constant 0, 0, [1, lo] ■ The entropy cost to dis- 
card the information obtained about each gas molecule 
appears to always exceed the entropy reduction achieved 
by the demon. Maxwell's demon is hence exorcised by 
Landauer's erasure principle. 

The second law stipulates that irreversible entropy 
production is positive in macroscopic systems. Ther- 
mal fluctuations are usually exceedingly small at these 
large scales and are therefore discarded. By contrast, 
fluctuations become predominant in microscopic systems 
and it has lately been recognized that the second law 
has to be generalized to properly take positive as well 
as negative entropy fluctuations into account This 
generalization takes the form of a fluctuation theorem, 
P(— S) = P{T,) exp(— S), for the probability distribution 
of the entropy production P(S) d, Q. Processes with 
negative entropy production can hence occur in small 
systems with small S, while they are exponentially sup- 
pressed in large systems with large E. The fluctuation 
theorem has recently been confirmed experimentally us- 
ing a colloidal particle in a modulated optical trap 
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impact of fluctuations on memory erasure in small sys- 
tems. By considering a generic system, we show that full 
erasure may be achieved in nanosystems with an entropy 
dissipation less than the Landauer limit. We express the 
Landauer bound in terms of the free energy difference of 
the erasure process and quantify the probability of hav- 
ing a dissipated entropy below that value. We moreover 
discuss an experimental set-up where our findings can 
be tested using currently available technology. To our 
knowledge, Landauer's principle has never been verified 
experimentally, despite its fundamental importance as a 
bridge between information theory and physics. One of 
the main difficulty is that in order to access the dissipated 
entropy, one needs to be able, like Maxwell's demon, to 
follow individual particles. Thanks to recent progress of 
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single molecule experiments, this is now possible 

The model. Following the original work of Landauer, 
we consider an overdamped Brownian particle in a one- 
dimensional double-well potential as a generic model of a 
one-bit memory Q . The state of the memory is assigned 
the value zero if the particle is in the left {x < 0) well 
and one in the right (x > 0) well. The memory is said to 
be erased when its state is reset to one (or alternatively 
zero) irrespective of its initial state. The potential barrier 
is assumed to be much larger than the thermal energy so 
that the memory can be considered stable in the absence 
of any perturbation. We describe the dynamics of the 
Brownian particle using the Langevin equation. 
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and a driven torsion pendulum 

Motivated by the recent development of nanotechno- 
logical memory devices 12, [l^, Ej], we investigate the 



where 7 is the friction coefficient and ^(i) a delta 
correlated Gaussian noise force with {^{t)) — and 
{^(t)i(.t)) = "2.0 5{t -t). The diffusion coefhcient is 
given hy D = ^kT. The double- well potential V{x) 
is of the standard form, V{x) — —ag(t)x'^/2 + 6x^/4, 
with a barrier height that can be controlled by the di- 
mensionless function g{t). The Brownian particle is 
additionally subjected to a driving force Af{t) with a 
driving amplitude A. In order to simplify the analy- 
sis of Eq. ([1]), it is convenient to introduce dimension- 
less variables defined as x = x/xm and t — at/^ where 
Xm = \/ a/b is the position of the (positive) minimum of 
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the potential [13] ■ The dimensionless potential is then 
V{x) = —g{t)x^/2 + a;**/4 and the rescaled driving am- 
plitude and diffusion coefficient are respectively given by 
A = A/aXm and D — D/^ax^. In the following, we will 
drop the bar signs and consider dimensionless quantities. 

As discussed in detail by Bennett in Ref. 0] , the mem- 
ory can be reset to one by i) lowering the barrier height 
and ii) applying an external tilt that brings the nanopar- 
ticle into the right well. We here characterize the erasure 
protocol with the help of the two functions g{t) and f{t) 
shown in Fig.[TJ The function g{t) defined as, 



1 - Csm{ijj{t ~ to)) if t e [to,to + tf] , 
1 otherwise. 



(2) 



lowers and raises the potential barrier in a time tf = t/2 
specified by the period, r = 2tt/u;, of the sine function. 
The parameter C controls the amplitude of the barrier 
lowering. On the other hand, the tilting function f{t) 
has a sawtooth shape parameterized as. 
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(i-io)/Ti ifte [tQ,tQ + Ti], 

l-{t-to- Ti)/t2 

if t e [to + Ti,h + tf] 

otherwise . 



(3) 



The two time constants ti and T2 verify tf = ti + T2 and 
are respectively the times during which the driving force 
Af is ramped up to its maximal value A and then down 
again to zero. The values of the parameters in Eqs. ^ 
and ([3]) are chosen in order to minimize the dissipated 
entropy. After a full erasure cycle, both g{t) and f{t) take 
back their original values and the double-well potential 
V{x) is restored to its initial shape. 

The energetics of the Brownian particle can be intro- 
duced by following the prescription of Ref. [lH . The vari- 
ation of the total potential energy is defined as AU = 
Uix{to + tf),to + tf) - Uix{to),to), where U{x,t) = 
V{x, t) — xAf{t) is the sum of the double-well potential 
V{x,t) and the driving potential —xAf{t). The work 
performed on the particle is in turn given by, 
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Since the driving force satisfies f{to + tf) ~ f{to) = 0, 
the work can be rewritten in terms of the functions g(t) 
and f(t) asW = //°+*^ dtx[Af(t) + {g(t) -l)x]. This ex- 
pression indicates tliat particles moving against the total 
force, F{x,t) = Af{t) + {g{t) — l)x, will generate work, 
while particles moving in the same direction will absorb 
work. This is the physical origin of work fluctuations in 
the present system. According to the first law, the heat 
dissipated into the bath is the difference between work 
and total energy change, Q — W~AU. It is important to 
note that the energy change during an erasure cycle van- 
ishes on average, (At/) =0. As a consequence, the mean 
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FIG. 1: (Color online) Time dependence of the two dimen- 
sionless functions g{t) and f{t), Eqs. ([2| and (|3}, that govern 
the erasure protocol. The function g{t) lowers the potential 
barrier, whereas f{t) induces a tilt that brings the Brownian 
particle into the right well (state one of the memory). 
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FIG. 2: Stochastic evolution of an ensemble of 100 trajectories 
during the erasure process. The state of the memory is zero 
or one with probability one half at initial time to- Erasure of 
the memory to state one is complete in case (a), while only 
partial in case (b). We have here used to = 0.01 and D = 0.02. 



work done on the system and the mean dissipated heat 
are equal, {W) = (Q). Moreover, the dissipated heat is 
related to the dissipated entropy via (Q) = TAS for a 
quasistatic transformation. The Landauer bound for the 
dissipated heat then follows as {Q) Landauer = fcTln2, or 
equivalently, {Q) Landauer = Dhi2 in dimensionless units. 

Numerical results. We have numerically integrated 
the overdamped Langevin equation ([T]) using the Heun 
method (Runge-Kutta n — 2) ^] with a Gaussian noise 
force generated by means of the Box-Muller method 21 1 . 
In our simulations, the particle is initially placed at po- 
sition a; = and is then left to equilibrate with the heat 
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FIG. 3: (Color online) (a) Memory erasure rate, defined as 
the relative number of particles ending in the right potential 
well, as a function of the amplitude of the tilting force Af{t); 
full memory erasure is attained for A > Ao — 0.15. (b) Mean 
work (W) and mean dissipated heat {Q) versus the tilting 
amplitude A for parameters uj — 0.01, C = 1.2, D = 0.02, 
To = 0.05tf and n = 0A5tf. 
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FIG. 4: (Color online) (a) Distribution of work and heat for 
a full erasure process, (b) Conditional heat distribution for 
particles initially in the left well. The existence of fluctuations 
below the Landauer limit is here clearly visible. Parameters 
are lo = 0.01, A = 0.2, C = 1.0(a), 1.2(b), and D = 0.02. The 
average work and heat is given by {W) = (Q) = 0.036 and 
the Landauer bound is {Q) Landauer — D\n2 = 0.014. 



bath at temperature T. After a thermalizing time to, the 
memory is therefore either in state zero or one with prob- 
ability one half. Two examples of an erasing sequence 
for different driving amphtudes (see below) are shown in 
Fig. m for an ensemble of 10^ particles. The erasure pro- 
tocol is applied at t = to for a duration of t/ . We observe 
that in the first case, all particles end up in the positive 
potential well. Full reset of the memory to state one is 
therefore achieved. On the contrary, only partial erasure 
is realized in the second case, as some Brownian particles 
remain in the negative potential well after the end of the 
erasure process. 

Figure [3] shows the erasure rate, defined as the rel- 
ative number of particles being in the right potential 
well at the end of the erasure protocol, as a function of 
the driving amplitude A. The simulations are performed 
with 10^ trajectories and a time slicing At — t * 5.10^^. 
We note that full erasure is obtained for driving ampli- 
tudes larger than a threshold value of ^o — 0.15. The 
mean work (W) and the mean dissipated heat (Q) dur- 
ing an erasure cycle are plotted in the inset: both in- 
crease monotonically with increasing driving amplitude. 
For the chosen parameters, the dimensionless diffusion 
coefficient is 13 — 0.02, corresponding to a Landauer 
limit of (Q) Landauer — -D lu 2 = 0.014. We stress that 
the mean dissipated heat always exceeds the Landauer 
bound, (Q) > {Q) Landauer, in accordaucc with Lan- 
dauer's principle. We have similarly studied the depen- 
dence of the erasure rate on the parameter C (not shown) . 

The full probability distributions of work and heat for 
a given protocol resulting in complete memory erasure 



are depicted in Fig. 21 Three important points are worth 
discussing. First, P{W) and P{Q) are broad distribu- 
tions with positive and negative fluctuations induced by 
the thermal bath. We note furthermore that the heat 
distribution is markedly broader than the work distribu- 
tion. The difference in width of the two distributions can 
be understood by noticing that, following the definition 
of Q, its variance is larger than the sum of the variances 
of work and energy, ctq > a^r + cr'^u > tr^. The sec- 
ond observation is that both distributions have a bimodal 
structure. The left peak arises from half the particles be- 
ing initially located in the right potential well, whereas 
the right peak stems from particles being moved from the 
left to the right potential well during the erasure process. 
The inset shows the conditional heat distribution pertain- 
ing to these particles. Third, and most importantly, we 
note that a significant fraction of trajectories lead to full 
memory erasure while dissipating less heat than the Lan- 
dauer limit, Q < (Q) Landauer- From thc conditional heat 
distribution shown in the inset, we find that 7.4% of all 
trajectories starting in the left potential well at time to 
yield a dissipated heat below (Q) Landauer- Full memory 
erasure in small systems can hence be obtained below the 
Landauer limit due to thermal fluctuations. 

Analytical considerations. We next derive the Lan- 
dauer bound for a nanoparticle in the double-well po- 
tential. Our starting point is the Jarzynski equality [l7| , 

{e-^^) = e~^^^ = ^ , (5) 

which relates the average exponentiated work done on 
the Brownian particle to the free energy difference AF = 
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—kTln Zi/Zq. Here Zq and Zi are the initial and fi- 
nal partition functions of the system. The Jarzynski 
equality is valid for arbitrary changes of the potential 
U{x,t) and has been successfully verified in biomolecule 
pulling experiments [23 • Before erasure, particles are 
thermally distributed over the entire double-well po- 
tential and the initial partition function is thus Zq = 
J^^dxexp{—PV{x)). After successfull erasure, all par- 
ticles are thermalized in the right potential well and 
Zi — dxexp{—l3V{x)) — Zo/2. As a result, the 
ratio of the partition functions is simply Zi/Zq = 1/2 
and AF — fcTln2. By using Jensen's inequality, 
(exp(— > cxp(— we immediately obtain, 
(W) = (Q) > fcTln2. The Landauer bound is thus 
given by the free energy difference of the erasure pro- 
cess and is only attained in the quasistatic limit when 
(W) = {Q) = AF. By further generalizing an argument 
presented in [2^ , the probability to observe a value of dis- 
sipated heat below the Landauer limit, {Q) Landauer — 1, 
with q > can be estimated to be. 



D = 10~^, as in our numerical simulations. The min- 
ima of the potential are then separated by Ax = 2xm = 
2^ ajh ~ 20 iim. These values are all com pat ible with 



Prob[(5 < {Q) Landauer " < exp(-/3g) 



(6) 



Small fluctuations below the Landauer limit are hence 
possible. However, large fluctuations, q ^ fcT, are expo- 
nentially suppressed, in agreement with the macroscopic 
formulation of Landauer's principle. We can therefore 
conclude that at the nanoscale, where fluctuations can- 
not simply be neglected, Landauer's principle has to be 
generalized in a way similar to the second law. 

Experimental implementation. The control of single 
nanoparticles in arbitrary two-dimensional force fields 
has been demonstrated in Refs. 15, In these exper- 



iments, fluorescence microscopy is combined with real- 
time feedback techniques to manipulate nanoscale ob- 
jects (from 50 to 200 nm) in solution via a position- 
dependent electrophoretic force. The force may be var- 
ied with high precision over nanometer distances and 
millisecond times. The trajectory of the particle is 
monitored with a high-sensitivity CCD camera. The 
Brownian motion of a 200 nm fluorescent polystyrene 
bead in a static double-well potential has been inves- 
tigated in Ref. [3] and the measured hoping rates be- 
tween the two wells have been successfully compared 
with Kramers theory. We propose to study memory 
erasure in nanosystems by extending the former exper- 
iment to a driven double-well potential according the 
erasure protocol discussed above. The work W done on 
the particle as well as the dissipated heat Q can be de- 
termined directly from the measured particle trajectory 
via Eq. By assuming a spacial diffusion coefficient 
Dx = 10/im^ s~^, a measuring time r = 60s and a tem- 
perature T — SOOif, we obtain the following realistic 
values for the parameters of the Langevin equation p]): 
7 = fcT/D^ ~ 4.10-1'' fcgs-i, a = -if It ~ ^.X^-"^ kg s''^ 
and h = Da^/kT ~ 40kgm~^s~^ with a dimension- 
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the experiments already performed in Refs. 

Conclusion. We have considered an overdamped Brow- 
nian particle in a double-well potential as a generic model 
of a one-bit memory. We have investigated the probabil- 
ity distribution of the work and the heat dissipated dur- 
ing an erasure process and demonstrated that full erasure 
may be reached by dissipating an amount of heat below 
the Landauer limit. We have shown that the occurrence 
of such events is exponentially suppressed, and there- 
fore not observable, in macroscopic systems. They play, 
however, an essential role in nanosystems and we have 
discussed a single-particle experiment where our predic- 
tions can be tested. Our main conclusion is that for small 
systems in general - and Maxwell's demon in particular - 
the macroscopic formulation of Landauer's principle does 
not hold, but has to be generalized to include heat fluc- 
tuations. 

This work was supported by the Emmy Noether Pro- 
gram of the DFG (Contract LU1382/1-1) and the cluster 
of excellence Nanosystems Initiative Munich (NIM) . 
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